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We show that the five-fold constraints due to (1) the observed nuclear modification of heavy 
quark jets measured via non-photonic electrons Raa (pt ~ 6 GeV) in central Au-I-Au collisions at 
200 AGeV, (2) the "perfect fluid" elliptic transverse flow of low transverse momenta pions, 1)2 (pt ~ 
1 GeV) reported for noncentral collisions, (3) the global pion rapidity density dN-^/dy, (4) the lattice 
QCD entropy density deficiency, S/Ssb, of strongly coupled Quark-Gluon Plasmas (sQGP), and 
(5) a causal requirement are analytically correlated in a class of gauge/string gravity dual models of 
sQGP dynamics. Current RHIC/BNL and lattice QCD data are found to be remarkably compatible 
with these models if the t'Hooft and Gauss-Bonnet coupling parameters lie in the range A ~ 10 — 25 
and < Xgb < 0.09. In addition, the observed five-fold correlation appears to favor color glass 
condensate over Glauber initial conditions within current systematic errors. 

PACS numbers: 25.75.-q, 11.25. Tq, 13.87.-a 



The combined observations of the quenching of hard 
(high transverse momentum or high quark mass jets) pro- 
cesses and the nearly "perfect fluid" eUiptic flow of soft 
(low momentum transfer) hadrons produced in Au+Au 
collisions at ^/s = 200 AGeV at the Relativistic Heavy 
Ion Collider (RHIC) l] have been interpreted as provid- 
ing evidence for the formation of a new form of strongly 
interacting quark-gluon plasma (sQGP) [3]. Further- 
more, bulk multiplicity (entropy) production systemat- 
ics have been interpreted as evidence for gluon satura- 
tion of the initial conditions as predicted by the Color 
Glass Condensate (CGC) model 0, However, it has 
been a challenge to find a single consistent theoretical 
framework that can explain simultaneously both soft and 
hard phenomena. These phenomena include 1) the nu- 
clear modification of high transverse momenta, pr > 5 
GeV, jet observables, 2) the bulk collective fiow observ- 
ables {pt < 1 GeV), and 3) the sQGP thermodynamic 
equation of state entropy deficiency relative to the ideal 
Stefan-Boltzmann limit as predicted by nonperturbative 
Lattice QCD (LQCD). Attempts to explain these based 
on weak coupling (perturbative QCD) parton transport 
approaches [a, S] , especially the surprising small vis- 
cosity needed to fit the elliptic flow data, require large 
coupling extrapolations as — gyj^j/AT: 0.5 — 0.6. 
At such large gauge couplings, on the other hand, the 
t'Hooft parameter A = gyM^c ^ 1 may already be large 
enough to validate string theory inspired AdS/CFT low 
energy approximations [3, . The great theoretical ad- 
vantage and appeal of gravity/gauge dual models 
is that they have led (for the first time) to analytic con- 
nections between a wide variety of thermodynamic and 
nonequilibrium dynamic variables at strong coupling that 
were not yet realized with traditional gauge theory tech- 
niques. 

In this Letter, we focus on predicted analytic con- 
nections between three fundamental properties of the 
sQGP: (1) its equation of state (entropy=S), (2) its 



long wavelength transport coefficients (viscosity^?]) , and 
(3) coupling between long wavelength near-equilibrium 
"soft medium" properties and short-wavelength non- 
equilibrated "hard probes" (energy loss per unit length, 
dE/dx). We show that these three properties together 
with global entropy and causality restrictions can provide 
valuable phenomenological constraints of higher dimen- 
sional gravity dual models of sQGP in heavy ion colli- 
sions. We consider here the constraints imposed by cur- 
rent RHIC and LQCD data on a class of gravity dual 
models that include quadratic as well as quartic cur- 
vature corrections to the classical Einstein-Hilbert ac- 
tion for the effective 5 dimensional dual gravity action. 
Therefore, we implement and extend the suggestion made 
in [l6[ to include perturbatively both the lowest order 
0{\gb 1/iVc) Gauss-Bonnet 7^2 and 0{l/\^/^) 
curvature corrections to the three properties above. 

Our analysis is based on the following remarkably sim- 
ple algebraic expressions relating these three fundamen- 
tal sQGP properties: 

S/SsB =i(l + AGB + ff|), (1) 

vis = ± (^1 - 4Agb + 15^) , (2) 

where s = S/V is the entropy density. The heavy 
quark jet relaxation rate, ^/tq, is controlled by /ig — 
^/XttT'^ /2Mq for a heavy quark with mass Mq in a 
plasma of temperature T . The relaxation time is related 
to the energy loss per unit length through tq(A, \gb) = 
— l/{d\ogp/dt) = -l/{dlogE/dx), where p — Mqjv 
and V = p/E. 

We note that the TZ^ correction 0{1JX'^^^) to the heavy 
quark jet energy loss is a new result [25| reported in this 
Letter and it is needed for a consistent additive perturba- 
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FIG. 1: (Color online) The gravity dual model correlation 
of the nuclear modification factor , Raaj oi hard (pr ~ 5.5 
GeV, 0-10% centrality) non-photonic electrons with the ellip- 
tic transverse flow moment, V2, of soft (pt = 1 GeV, 20-60% 
centrality) pions assuming Glauber (left) and CGC (right) 
initial participant distributions and dA'^/dy = 1000 are com- 
pared. The green region corresponds to the range of entropy 
function deficiencies, S/Ssb, consistent with infinite volume 
extrapolations of lattice QCD data [sj. The red (blue) el- 
hpse indicates PHENIX (STAR) 0-10% centrality R'Xa and 
20-60%U2 centrality Au-^Au 200 AGeV data from RHIC 
[s^ l- The dashed blue curves show fixed S/Ssb = 0.8 while 
for the solid black and dotted-dashed curves S/Ssb = 0.85 
and S/SsB ~ 0.9, respectively in the range A = 5, ...,30. 
The causality bound curve with Xgb = 0.09 is indicated by 
solid red. The gray ellipse shows preliminary V2 for minimum 
bias data with estimated non-flow effects subtracted 13911. 



tive application of this class of AdS/CFT models to heavy 
ion reactions. We consider here only heavy quark observ- 
ables because gravity dual string drag models apply 
only for heavy Mq > ^fXT jets. The heavy quark jet 
drag is then modeled as a trailing string moving in a black 
brane background according to the classical Nambu-Goto 
action (with dilaton neglected) Ang = "^TrtF / '^^'^V^ 
where g = det gab = G^j.vdaX^^dbX" is the induced 
worldsheet metric, cr" ~ (r, a) are the internal world- 
sheet coordinates, G^^{X) is the background metric, and 
= X^^ (r, a) is the embedding of the string in space- 
time. The trailing string ansatz (where t = t, a = u 
and X^^{t^u) ~ {t,xo + vt + £^{u),0,0,u)) describes the 
asymptotic behavior of a string attached to a moving 
heavy quark (the string endpoint) with velocity v in the 
X direction and located at a fixed AdS radial coordinate 
Um 3> Uh [13*1 . The black brane horizon coordinate Uh oc 
Ta' is is determined by G'oo('W/i) = 0. Using the ansatz 
above and the string's classical equations of motion, one 
can show that the drag force dp/dt = —Cvj (27ra'), where 
C is a constant determined by the negativity condition 



g(u) = Gun (Goo + v^G^^) ( 1 + 



GmGrr 



<0 (4) 



for < u < Um- However, both the numerator and 
denominator in Eq. (0]) change their sign simultaneously 
at a certain u* [13] given by the root of the equation 
Gooiu*) + «2c(u*) = 0. This fixes C = G,.(m*) 
and dp/dt — ~v Gxx{u*)/ {2TTa'). Neglecting higher- 
order derivative corrections in J\f — 4 SYM one finds 
u* = Ufiy/j, where 7 = l/\/l — v^. The condition that 
u* < Um leads to a maximum "speed limit" for the heavy 
quark jet to be consistent with this trailing string ansatz 
given by 7n 



< u'i/ul 



14| 



Using the metric derived in [1C| to 0{a ), one can 
compute the effects of quartic corrections on the drag 
force [3l and determine u* perturbatively to 0(A~'^/^) 



as [26| : 
drag force 
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Thus, the 



The heavy quark mass at T = is AIq — 

and, to leading order in 1/A, u^/u\ ~ 
corrected u* displayed above defines a new speed limit 

7m ~ TE (^f^ ~ ' after neglecting terms of 

0(1/7, 1/iVc). Note that 7™ and dp/dt decrease with 
increasing 77/ s. 

For our applications, we consider the range A ~ 5 — 30 
and |Agb| < 0.1 with fixed — i. In this parame- 
ter range the 1/A"^/^ and Xgb ^ l/-^c corrections are 
comparable. We neglect known (but formally) higher or- 
der terms [l^l 0{y/\/Nl) in this first attempt to test 
predicted dynamical correlations between hard and soft 
phenomena in high energy A-fA collisions. 

The small Gauss-Bonnet parameter Acs = (c — a)/4c 
is related to the central charges c and a (related to the 
conformal anomaly in curved spacetime) of the dual CFT 
as noted in Eq. (2.14) of Ref. [l6|. Varying Xqb provides 
a parametric way to explore deformations of the original 
A/" = 4 SU{Nc) SYM theory. Interest in Gauss-Bonnet 
deformations were heightened when Kats and Petrov 21 1 
argued that for TV = 2 Sp{Nc), Xgb = l/8Nc, the KSS 
viscosity bound on r]/s > l/An was violated by 17% 
for Nc = 3. As further shown in [l^, a large class 
of other effective CFTs are now known to lead to sim- 
ilar Xgb oc 1/A''c effects. However, the analysis of Refs. 
[ll]-[20| revealed that Xgb deformations are limited by 
requirements of causality and positive energy flow to a 
narrow range —7/36 < Xgb < 9/100. 

In order to convert tq into the observed nuclear mod- 
ification of single non-photonic electrons, i?^^(pT — 
5.5 GeV), from quenched heavy quark jets, we follow 
and extend pTj by using here the generalized drag force 
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FIG. 2: (Color online) Five fold phenomenological constraints 
in the t'Hooft and Gauss-Bonnet parameter space {X,Xgb)- 
(1) The green region with black dashed contours is from lattice 
QCD constraints on 0.8 < S/Ssb < 0.9. (2) The cyan region 
is determined from noncentral elliptic flow V2 (pr = 1 , 20 — 
60%) — 0.11 ± 0.01 38]. Blue dashed contours correspond 
to fixed Anrj/s = 1,1.5,2. The inversion Agb(A, y2(v/s)) is 
based on minimum bias viscous hydro results of [3l| assuming 
CGC initial eccentricities scaled by a factor 1.1 in panel (a) 
and unsealed 1.0 in panel (b). The entropy is constrained 
by the (3) dA'^/dy — 1000 pion rapidity density in central 
collisions. (4) The gray region and contours are determined 
from central Raa{Pt = 5.5 GeV) = 0.25 ± 0.07 data [H. 
(5) The horizontal red line constraint is the causality upper 
bound for Ags"= = 0.09 [li,!!^]. The yeUow trapezoidal region 
with the purple boundary is the intersection of the five fold 
constraint bands. Note the red circled five fold conjunction 
area in panel (a) (A ~ 13, \gb ~ 0.08) that is absent in the 
unsealed panel (b). 



in Eq. ^ to compute the path length, L, dependent 
heavy quark fractional energy loss e(L). The heavy 
quark jet nuclear modification factor is then Raa — 



((1- 

tral index 



!l, where nQ{pT) is the flavor dependent spec- 

— I obtained from 
I as used in 
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FONLL production cross sections [2: 
The path length average of the nuclear modification at 
impact parameter b is computed using a Woods-Saxon 
nuclear density profile with Glauber profiles Ta{x±) with 
ctnn — 42 mb. For 0-10% centrality triggered data 
both Glauber and CGC geometries lead to similar nu- 
merical results [iSl- The distribution of initial hard jet 
production points at a given x± and azimuthal direc- 
tion (j) is taken to be proportional to the binary par- 
ton collision density, TAA{x±,b). We assume a longi- 



tudinally expanding local (participant) parton density 
pix±,b) = XPpartix±,b)/T, whcic X = idN^/dy)/Np^^t 
and Ppart is the Glauber participant nucleon profile den- 
sity. However, we evaluate ([3]) with a reduced tempera- 
ture TcFT = 0.1A{S/ SsbY^^Tqcd to take into account 
the smaller number of degrees of freedom in a strongly- 
coupled QCD plasma, which is similar to the prescription 
given in 30] . We compute the heavy quark modification 
factor R^A ^'^^ 



RiAiPT,b) = J\cl> J d'x^^ 



TA Aix±,b) 

iVbin{6) 



X exp 



dT 



Tc{x±_ + Te{(l)), 



(6) 



where iVbin is the number of binary collisions. Here, 
To = 1 fm/c is the assumed plasma equilibration time 
and Tf is determined from T{l,Tf ) = Tf = 140 MeV, 
i.e, the time at which the local temperature falls below 
a freeze-out temperature taken from ^Sl*] . Systematic er- 
rors associated with the freeze-out condition will be dis- 
cussed elsewhere. 

In order to compute W2(l GeV) for the C(20-60%) cen- 
trality class, we employ a linear fit to the numerical vis- 
cous hydrodynamic results of Luzum and Romatschke 
31 ] . The dependence of V2 on viscosity for both Glauber 
32] and CGC 4.] initial transverse profiles can be well fit 
with 



V2{PT, ri/s,C) = a{pT) e2(C) (1 -brj/s) 



(7) 



where e2(C) = (y^ — x^)c/{x^ + y^)c is the average initial 
elliptic geometric eccentricity for the centrality class C. 
To rescale the minimum bias viscous hydro results of Ref. 
3l\ to the considered 20-60% centrality class we use the 



factor e^'""''(20-60%)/e^''"'^(0-92%) = 0.317/0.281 
1.128 from Ref. [s^]- Our fit to the rescaled numerical 
results of [lH give b « 2.5 and a{pT = l)e2(20 - 60%) « 
0.14 (0.098) for CGC (Glauber) initial conditions. We 
consider this 20-60% centrality class because, as shown 
in Fig. 23 of , there is good agreement at pt '^-^ 1 GeV 
between STAR W4) and PHENIX V2{BBC) data and 
non-flow effects [39| that complicate the interpretation 
of minimum bias data in (sij are reduced. 

We find (see Figs. [T] and [J) that a (reasonable) com- 
bination of model parameters (A and Xgb) can account 
for the correlation between the reported R\a 0; S| 
and V2 [3^ taking into account the LQCD constraint 
on the equations of state deficiency S/Ssb [13] (in the 
range T ~ 2 — 4Tc). Our results suggest that within 
current systematic errors a small positive quadratic cur- 
vature correction with < Acs < 0.1 is preferred phe- 
nomenologically. Moreover, CGC initial eccentricities are 
favored over Glauber eccentricities for the centrality class 
considered. However, as emphasized by Fig. 2b, the re- 
sults are rather sensitive to the Ags(A, W2('7/s)) inver- 
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sion. We find that a 10% reduction of tiie viscous hydro- 
dynamic V2{ri/ s) (as for minimum bias centraUty) virtu- 
ally eliminates the yellow overlap region and could falsify 
the AdS/CFT description based on Eqs. (1-3). Great 
care is called for at this time to avoid premature con- 
clusions. Systematic studies on viscous hydrodynamic 
dependence on centrality cuts and improved experimen- 
tal control over non-flow corrections will be needed be- 
fore definitive conclusions could be reached. The good 
news demonstrated by comparing Fig. 2a and 2b is that 
if improved theoretical and experimental control (bet- 
ter than 10%) over elliptic flow systematics and the nu- 
clear modification of heavy quark jet observables can be 
reached, then rather strong experimental constraints on 
the AdS/CFT gravity dual model parameters could be 
achieved. We close by emphasizing [27| that future com- 
parison of the nuclear modification of identified bottom 
and charm quark jets at RHIC and LHC combined with 
the fivefold (hard/soft) constraints considered in this Let- 
ter will provide especially stringent tests of AdS/CFT 
gravity dual phenomenology applied to high energy heavy 
ion reactions. 
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